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Abstract 

One-dimensional quantum scattering from a local potential barrier is considered. 
Analytical properties of the scattering amplitudes have been investigated by means 
of the integral equations equivalent to the Schrodinger equations. The transition and 
reflection amplitudes are expressed in terms of two complex functions of the incident 
energy, which are similar to the Jost function in the partial-wave scattering. These 
functions are entire for finite-range potentials and meromorphic for exponentially 
decreasing potentials. The analytical properties result from locality of the potential 
in the wave equation and represent the effect of causality in time dependence of the 
scattering process. 

PACS numbers: 0230, 0380, 1155, 0365. 

Short title: Analytical properties of scattering amplitudes 
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1 Introduction 



The problem of the tunneling time has been attracting a considerable attention for decades 
[jl], 1^, 1^, ^. It is indeed important to understand the effect of causality upon the particle 
and wave propagation. The problem is that, in the conventional time-independent for- 
malism, the causality manifests itself indirectly, namely, in analytical properties of the 
transition amplitudes as functions of (complex) energy. The relations between causality 
and analyticity were under an intensive investigation in the 60-ties, when the concept of 
the S'-matrix dominated in particle physics[^, ^. At that time, however, the analysis was 
aimed mainly at three-dimensional scattering processes, for central-symmetric potentials 
in particularly, |^. Partial-wave scattering amplitudes in the complex energy plane were 
considered also in the theory of multi-channel nuclear reactions 0. Analytical properties 
of the scattering matrix have been used also in other fields. For example, in the theory of 
multi-terminal mesoscopic conductance the properties of the multi-probe S'-matrix were 
used[|lO| to obtain the conductance and to establish the time ordering of the incoming and 
outgoing lead states. In another work|TT[|, the low-frequency behaviour of dynamic con- 
ductance was related to the phase-delay times for the carrier transmission and reflection, 
which are given by the energy derivatives of the S'-matrix elements. The most imporatnt 
feature of the transition amplitudes for various physical processes in the energy represen- 
tation is their analyticity in the upper half of the complex plane. Sometimes one can find 
out more about singularities in the lower half-plane, investigating dynamical equations 
specific for the physical problem, like the Schrodinger equation for scattering. 

The one- dimensional 'scattering', i.e. the potential-barrier problem, is somewhat more 
complicated than the scattering off a force center, since the system has two channels, 
corresponding to two waves running in the opposite directions from the potential region 
in the final state. Thus instead of one analytical function, the partial wave scattering 
amplitude, one deals with two analytical functions, transmission and reflection amplitudes. 
(The unitarity condition holds in both cases.) The familiar arguments of the scattering 
theory must be extended properly to the one-dimensional case. The analytical properties 



of the one- dimensional S matrix were discussed, in particular, by Faddeev|T^, [T^ and 
Newton fl^ in view of the inverse problem. The singularities in the complex energy plane, 
caused by bound and quasi-bound states were considered more recently [16|, [l^ . 



The purpose of the present work is to investigate the location and character of singu- 
larities in the scattering amplitudes, owing to the potential shape. The investigation is 
based upon the Schrodinger equation with a local potential. That enables one to reveal 
general features of the amplitudes, depending on the character of vanishing of the poten- 
tial outside its domain. The analytical properties are essential for the application to a 
space-time picture of the barrier transmission. 

In Section 2, the 2x2 scattering and transition matrices are introduced and related to 
the resolvent of the Schrodinger operator. Next, in Section 3, two complex functions are 
defined for the potential barrier problem, which are related to elements of the monodromy 
matrix. Their role is similar to that of the Jost function in the S'-wave potential scattering. 
These functions have nice analytical properties, which are proven in Section 4 by means 
of the Volterra-type integral equations. It is shown, in particular, that the functions 
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have no singularities in the whole complex energy plane (except for the infinity), if the 
potential has a finite range. Singularities appear if the potential behaves exponentially in 
the asymptotics, and the slope of the exponent determines the distance to the singularies 
nearest to the real energy axis. Some examples are given in Appendix. 



2 Transition operator and the /S-matrix 

The evolution operator is given by the Laplace transform of the resolvent of the Hamil- 
tonian if, 

exp{-ttH) = f G,e-'''de, = {H - s)-\ (1) 

Here Too is the contour in the complex e-plane, running from — oo to +00 above the real 
axis, where the singularities are situated. In scattering problems, the Hamiltonian is a 
sum of a potential operator and the kinetic energy, having the eigen-vectors describing 
free particle states, 

H = Ho + V, ^o|k>=e(A;)|k> . (2) 

In the nonrelativistic scattering theory, we shall use the units where e(A;) = k"^. The 
transition operator is introduced as follows, 

G', = G'f -GftGf , G^^^ = {Ho-e)-\ (3) 

It can be expressed directly in terms of the resolvent G^, 

t = V- VG,V, (4) 

As follows from the time- inversion symmetry (the reciprocity principle), 

<k|T;|ko>=<-ko|f,|-k>. (5) 

It is easy to see, using the standard definition of the scattering operator S, that its matrix 
elements are expressed in terms of the transition operator on the energy shell, 

< k|,S|ko >= lim < k\ei^^°e-'*^ei*"°\ko > 

t— »oo 

- 5{k - ko) - 27riS[e{k) - e{ko)]T,,,{e) . (6) 
Here the scattering amplitude is given by 

T,,„(£) klf.Iko >, e{k) =e = e{ko), (7) 

where u = k/k, and the standard normalization is used: < k|ko >= 6{k — ko). 

In the one-dimensional case k = uk, where i/ = ±1, corresponding to two possible 
directions of motion for a given energy, and 

S{k — ko) = vS[e{k) — e{ko)]Si,^fj, v — de/dk. (8) 
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The elements of S and T (on the energy shell) are given by 2 x 2 matrices S and T, 

27r? 

<k\S\ko>=v6[e{k)-e{ko)]S,,,, S = I - —T, (9) 

By definition, if S exists, it is a unitary operator. This fact implies a unitarity condition 
on the scattering amplitude, which reads 

SS^ = I, —{T-T^) = --TTl (10) 

As we will show, the analytical properties of T{e) follow from locality of the potential, 
and Eq. (4). 



3 General properties of the transition amplitudes 

We consider the Schrodinger equation, 

Hip = K^^, H = -{d/dxf + V{x), (11) 

where — cxo < a; < cxd, and the potential is local, i.e. |V"(x)| = o{l/x) as \x\ —>■ oo. In the 
coordinate representation, the resolvent can be expressed in terms of two fundamental 
solutions of the Schrodinger equation, y±{x), satisfying the proper boundary conditions 
at ±oo, respectively. 



< x|Ge|xo >= ^ — , a;>/< = max/mm(x, Xo), (12 



y -{x<)y+{x> ) 

y±{x) ^ 6=^*"'^, as X -> ±oo, (13) 
w{y^, y+) = y'_y+ - y-y+ = const. 



As soon as £ = is introduced in the Laplace transform (1) for the upper half of the 
complex plane, the solutions y± defined above vanish at ±oo, respectively. Thus one 
has the properly defined resolvent for the elliptic operator H satisfying the Sommerfeld 
radiation condition at infinity ||18|. 

For large where the potential vanishes, the asymptotics of the fundamental solu- 
tions are given by 

y_(x) = ae"*'''' + fee*"'' for x -> +oo (14) 
y+{x) = fo'e"*"'' + ce*"'' for x -oo. 

In principle, the solutions are defined by (13) for Ken > and Imn +0, yet the 
analytical continuation to the whole complex K-plane is considered in the following. For 
real potentials V{x), the complex conjugate functions y±{x) are also solutions of the 
Schrodinger equation, satisfying the boundary conditions conjugate to (13). Calculating 
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the Wronskians (which are independent of x) at x ±00 for various pairs of the solutions, 
one gets a number of relations between the complex parameters a, b, b', c: 

w{y^, ?/+) = — 2ma = — 2mc ^ a = c, (15) 

w{y^,y^) = const ^ |ap — |6p = 1, (16) 
w{y+,y+) = const ^ |c|^ — = 1, 

w{y^,y+) = const b' = —b. (17) 

Expressing y± in terms of the fundamental solutions, one has 

b I lb 

y- = -y^ + -y+, y+ = -y- - -y+- (18) 

a a a a 



The analytical continuation of these solutions to the complex k plane, by y±{x) = y^'^{x) 
implies a symmetry of «(«) and 6(k) with respect to the imaginary /t-axis. 



a{K) = a{-R), b{K)=b{-R). (19) 

Thus, the asymptotics of the solutions depend only on two complex functions a(K) and 
6(k), satisfying one real condition (16), and subject to the symmetry (19). The transition 



amplitudes, elements of the S matrix and T matrix, and of the monodromy matrix |T9 
are given in terms of these two functions. 

If the potential is displaced, b gets a phase shift, 

V{x) ^V{x-d) -^a-^a, b^ be''^'^'^. (20) 

For symmetric potentials one gets an additional relation, 

V{x) = V{—x) ^ y_(— x) = y+(x) : a = c, b' = b, (21) 

so that b is pure imaginary, in view of (17). 

As soon as the resolvent is known from (12), the elements of the transition operator 
in the momentum representation are obtained immediately, by (4), 

< A;|f,|A;o >= V{q) jj dxdxoe-'^''^'^°''°V {x) < x|G,|xo > V^(xo), (22) 
where q = k — k^, and 

1 f+00 

V{q) = — / Vix)e-"'''dx. (23) 



27r 

The double Fourier transform in (22) is performed by means of the Schrodinger equation 
Vy = y" + K^y, leading to the following integrals 



-ikx 



+ (^2_^2) r e-'^VlO^e, (24) 

J —00 

e-^'^ViOy+iOd^ = -« + zA;77+)e-^^^ 

H^'-k') e-^'^vAOdt (25) 



where we have introduced the functions r]±{x) vanishing at ±00, 

r]^{x)=yi{x)-e^'^^. (26) 

The result is 

< k\t\ko >= r dxe-'^^ix) \iK - l)y4x)e-'^^ + («: + hy4^)e'' 

iTXiw J-00 L 2 2 
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/oo 
-00 



-{k^ - A;2)(«:2 - A;2)— yy dxdxoe-''^+'''°^°g{x,Xo), (27) 

where w = — 2ma, and 

5f(x,a;o) = ?7_(a;<)r7+(a;>)/u;. (28) 

On the energy shell where fc^ = = fcg, the contributions from the integrals vanish 
in Eqs. (24-27), and one has 



2™ \ (3 a r a\ -h I 



where 



a= dxe-'^''V{x)y+{x) = dxe"'''V{x)y4x), (30) 

J — 00 J — 00 

/OO roo 
dxe'''''Vix)y+{x) = / dxe-'^'^V{x)y-ix), 

a = l-—, b = ^. (31) 
2m 2m ^ ^ 

The functions «(«;) and /5(k) are free of a pole at = 0, and are related by the unitarity 
condition, 

a — a = —{aa — [3[3). (32) 

The transmission and reflection amplitudes are 5*++ = 1/a and S'+_ = 6/a, respectively, 
and det S = a/a. (The latter equality supports the analogy of a to the Jost function. If 
there is no reflection, 6 = 0, then a = e"*'', where 6{k) is a real phase shift.) Besides, if 
the potential is even, the matrices S and T are symmetrical. 

Note that because of the relations (15)-(17) the monodromy matrix composed of 
a, b is quasi-unitary, 

MEM^ = E, (33) 
Under a displacement of the potential, Eq. (20), M is transformed to 

M ^ DMD\ D = [ J ) • (34) 
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If V{x) = Vi{x — di) + V2{x — ^2), d2 > di, and the potential domain consists of two 
intervals separated by a forceless gap, one has the following superposition rule, 

M = D1M1DID2M2DI : 
a = aia2 + hih2e^'<'^'-'^'\ b = ai^se"'*'^'^^ + d2he-^'^'^\ (35) 

which is a consequence of Eqs. (14) and (20). 

4 Integral equations and the analytical properties 
4.1 Volterra equation and series solutions 

Analytical properties of the transition amplitudes can be derived from the integral equa- 
tion satisfied by the fundamental solutions y±{x), 

y±{x) = e±-- + - r sin«(x - 0^(Oy±(0^e- (36) 

K J±oo 

These equations arc of the Volterra type, so the solution exists and admits an analytical 
continuation to complex k for local potentials. In order to separate asymptotic oscillations 
of the solutions, let us introduce new functions (the integrals are evaluated by Eqs. (24-25) 
for k — ±«;), 

A_{x) = r e'-^ViOy-iOd^ = {y- - ^«y-)e^'^" + 2iK, 
J —00 

B-{x) ^ r e-'^^ViOy-iOdC = {y- + i^^y-)e-'^', (37) 

J— 00 

A+{x) = / e-'^^ViOy+iOdC = -{y+ + ^/^t/+)e-''^^ + 2ik, 

Jx 

roo 

B+{x)^ e^^^v{Oy+{OdC = -{y+-^^y+)e'^''- (38) 

Jx 

It is easy to see that 

y±{x) = (1 - ^M^)) e*'" + 5^^±(^)e'''". (39) 

so A and B have the definite hmits at x — > ±00, cf. Eqs. (30), 

a = A+{-oo)=A_{+oo), l3 = B_{+oo) = fi+(-oo), (40) 

A±(±oo) = = 5±(±oo). 

The pairs of functions {A, B) satisfy a system of first-order differential equations with 
zero initial conditions at infinity. Setting the equations into the integral form, one gets 
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from (37), in particular, for A- and 

A J.) = £ V(() (l - 1-A4() + ^S-tt).-0 '^^ 
B-C.) = £ V(f ) (d - 5^^-4-(f + 5i^B-(f )) (42) 

This form is especially suitable for the perturbative expansion of a and (3, namely, 

/ + 00 
-oo 

+ - r^di2 t'di,V{i2)V{ii)e'^^^^-^-^ sin «:(6 - 6) + " " " , (43) 

/t^ J —oo J —oo 

/ + 00 
r(Oe-^^^«rfe 
-oo 

+ - r^d^2 /^'^^ei^(6)m)e-^''^^^+«^^sin«:(6-ei) + ---- (44) 

K J —oo J —oo 

Using these expansions, one gets a sort of the Pade approximation for the S'-matrix in 
(29). Note that as k ^ 0, one has a — (3 —>■ while a and f5 are regular, so expanding 
in powers of k one can get, for short-range potentials, an analogue of the effective-range 
approximation . 



4.2 Analytical properties 

The perturbative series are also used to prove the analytical properties. First of all, one 
can extend to Eq. (36) the standard arguments of the scattering theory which are 
based upon the inequality 



sin /€(x - e) I < C '""7 ^ \im.\{x-o^ (45) 









1 + 


\kx\ 



where x > C,, and C is a constant. Thus one proves that y±{n) are analytical in the 
domain in the complex /t-plane where 



f exp[-^{\lmK\ ±lmK)]V{^)d^ < oo. 

J —oo 



In particular, if V{x) = for x < x_, for some x_, there is no irregularity as x — — oo. 
Similarly, the limit x — > +oo is considered. The fundamental solutions are analytical in 
K, as soon as these two limits are regular. 

One may rather modify the method and apply it directly to the functions we are 
interested in, given by Eqs. (41-42). The substitution 

A4x) = f{x)e'^^^\ B4x) = g{x)e-'^^^\ w{x) ^ ^ T d^ViO, (46) 

2k J-oo 
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eliminates the diagonal terms in the differential equations, and they are reduced to 

f'^fo + P+{x)g, /o = y(x)e-^-W, (47) 

g' ^ go + P_(x)f, go^V(x)exp[-2iKX + iw(x)]. (48) 

Here 

and the initial conditions are /(— oo) = = g{—oo). Note that f{x) and g{x) have limits 
as a; ^ oo, which are a and f3, up to conjugate phase shifts, provided that the potential 
is integrable, and w{oo) is finite. The solution to Eqs. (47-48) is given by the series, 

/=f:/n, fi{x)=r dMo, Ui{x) = r d^P409n{o, m 

jj^—l J —oo J —oo 

g-Egn: gi{x) ^ r d^goiO: gn+l{x)^ r d^P4^)fn{0- (51) 

n=l 

Upper bounds for and gn for positive and integrable potentials, can be obtained by 
iteration. Prom Eqs. (50-51), one gets for even n. 



n— 1 



P+(en-l)P-(en-2)...P+(ei)^o(eo) I] 

m=0 

5n(x) = / P_(en-l)P+(en-2)...P-(ei)/o(eo) R ^U, (52) 

m=0 

and for odd n, 

Ux) = f p+(e„-i)p-(en-2)...p-(6)/o(eo) n ^u, 

m=0 

gn{x) = / P.{^n-l)P+{^n-2)...P+{^l)go{^o) U dU- (53) 



n-1 



m=0 



The integrations take place in the domain A„ : — oo < .^o < < ■ ■ ■ < Cn-i < x. If we 
assume that V{x) > 0, kw{x) is a real, bounded and non-decreasing function of x. For 
Im K > 0, we shall use the following inequalities, 

\P4^i)gom < ^l/o(6)l, 

|/o(^o)| < ^(^o), 

\P-{0\ < ^|e-2«€||e2™(-)|, ^<x. (54) 
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Similarly, for Im k < 0, 



|P-(6)/o(6)| < ^l^o(^o)|, 

\P+iO\ < ^|e2^'^«||e-'^"'(^)|, C<x. (55) 
Using these inequalities one can show that, in the upper half of the complex K-plane, 

l/„Wi < |2.|M£)£, 

n! 

b„(x)e-^^-(^)| < |2«:|M^|^_(^)|, (56) 
while in the lower half of the complex K-plane, 

|/n(x)e2-(-)| < \2k\^-^^^^\u4x)u4x)\ 

\9n{x)\ < \2k\LSJ1-^\u_{x)\, (57) 

u±{x) = ^ r rfeV^(0exp(±2^O. (58) 

Ik, j-oo 

It is assumed that the integrals exist for real k and have definite limits as x — > oo (the 
Fourier transform of V). 

For positive and integrable potentials, f{x) and g{x), and thus «(«;) and /?(«;), are 
given in terms of infinite series which are absolutely convergent in the domain where 
the corresponding Fourier transforms of V{x) exist, Eq. (58). Thus a{n) and a{K) are 
analytical in the upper half K,-plane. The singularities of q;(k) and /9(k) appear if the 
integrals in Eq. (58) diverge as a; — > oo, which may happen only at finite distances below 
the real axis. 



where 



4.3 Finite-range potentials 

The functions a{n) and /?(«;) are entire, if V{x) = outside an interval (a;_,a;+). The 
analyticity is an immediate result of Eqs. (56-58), as u± are finite for potentials with a 
finite support. 

A more direct proof, as well as additional information on the analytic continuation 
into the complex re-plane, can be obtained from explicit expressions for a(re) and 6(re). Let 
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us introduce two real (for real k) solutions of the Sclirodinger equation, zo{x) and zi{x), 
specified by the following initial conditions, 

Zoix_) = 1, Zq{x^) = 0, 

zi(x_) = 0, z[{x-) = l. (59) 

From the continuity of the wave function and of its first derivative ai x = x±, one gets 
the following expressions for a and b, 

„ift(x_|_— x_) 

[-Co + ^^i(o + C[) + ^%], 
[Co-HCo-C[) + f^\i], (60) 



2iK 



2iK 

where Co,i = ^o,i(2;+)- As soon as ( and (' are analytical functions of k^, by the Poincare 
theorem[0, q;(k) and /3(/t) are also analytical in the whole complex K-plane. 
For large \k\ and smooth V{x), one can use the semi-classical approximation, 

I p_ 

Co = W — cos6', Co = -VP-P+sin6', 
V P+ 

Ci = , ^ sing. ([=.^cose, (61) 

V 

where 

9{k) = j - V{x)dx, p± = ^k'^ - V{x^). (62) 

In this approximation one gets 



6 

a = [(k^ + P-P+) sin 9 + zfi;(p_ + p+) cos^^l, 

2z«;v^IpT 

6 = [(k^ — P-P+) sin 6* — ift(p_ — p+) cos 6]. (63) 

2zK^/pIp7 

Note that this result is exact for the square- well barrier, Eq. (68), and the analytical 
continuation to the complex plane is possible. Asymptotical locations of zeroes of a in 
the complex K-plane are given by the equation 

exp|-2.0(.)| = '."-^-^l^-'i (64) 
Evidently, there are no zeroes in the upper half-plane for V{x) > 0. 
4.4 Exponentially decreasing potentials 

If V{x) oc exp(=l=2s±x) as x ^ ±oo, {s± > are constant); «(«;) and /?(k) are no longer 
entire functions. The singularities appear, when V{x) is not small enough to suppress 
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exp(±2mx), so the integral in Eq. (58) is diverging as x — >■ ±00. The singularities nearest 
to the real axis appear at Im k = —s± for / and and at Im k = ±s± for g and 

Explicit expressions for a(/t) and revealing their singularities, can be obtained, 
assuming that 



V{x) = e2«-(^-^-), x<X. 

V{x) 



^2 -2s+{x-x+)^ a;>a;+. 



(65) 



where v± are constants. The solution to the Schrodinger equation for a;_ < x < x+ is 
still a linear combination of zq and zi, while for x < X- and x > x+ it is given by linear 
combinations of the appropriate Bessel functions. Using the matching conditions for the 
wave function at x±, one gets. 



2m 



-r(l + z/+)r(l + z/_)(a+/2)-'^+(a_/2)- 



2m 



-r(l + z/+)r(l-z/_)(a+/2)-'^+(a_/2)^ 



CoJu+{cT+)J-u^{cr_) + zi;+CoA+(o-+)^-i._(o-_) 
+iv^CvJu+{(y+)J' - v+v^CiJ' v+{(y+)J' -u-{(^-) 



(66) 



where z/± = —iK/s± and (t± = iv±/s±. The singularities are due to the T functions, as 
soon as {a /2)^^ Jy^a) is known pO| to be an entire function in both a and z/, while C, and 
are analjd;ical functions of k^, by the Poincare theorem. Thus,both a{K) and h{K) have 
infinite series of equidistant poles on the imaginary axis: at k = —ins± for a{K) (the poles 
are double if s_ = s+), and at k = ^ins± for 6(k) (where n is any positive integer). 

The minimal distance of the singularities from the real K-axis, that was derived from 
the integral equations, is non-zero for every potential with an asymptotic exponential 
decline. The results of Eq. (66) are less general. If the assumption of Eq. (65) is relaxed, 
the poles can move off the imaginary axis, as one can see in Eq. (72) below. 



4.5 Singularities of the S matrix 

The singularities of the T and 5" matrices are of physical importance, when the time 
dependent process is considered. They are given by zeroes of a{K) as well as by the poles 
of &(k), which do not coincide with those of a{K). 

As was proven in section 4.2, for positive and integrable potentials, a{K) is analytical 
in the upper half K-plane, which is a result of causality; (3{k), and thus &(k), may have 
singularities at finite distances above and below the real /t-axis. 

The pattern of singularities of a{K) and 6(k) in the complex K-plane is determined by 
the asymptotic decline of the potential. For potentials, having an asymptotic decline faster 
than exponential, (e.g. finite range potentials, and the Gaussian barrier), no singularities 



12 



appear for finite k. For potentials with exponential asymptotics, the distance from the 
real K-axis to the nearest singularities is determined by the slope of the potential at ±00. 

It is important that a{K) has no zeroes in the upper half plane. As soon as the function 
is analytical, the number of its zeroes is given by the integral 

1 [da 
N = — 6 — , 
27rz Jc a 

where the contour C encloses the upper half of the complex plane. For positive and 
integrable potentials a^n) is limited in the upper half plane, so a^n) — > 1 as k ^ 00, and 
the integral is zero. 

The location of zeroes of a{K) in the lower half of the complex K-plane, depends on 
the specific barrier considered. For finite-range potentials the asymptotic distribution of 
zeroes is given by Eq. (64). Other examples are considered in the appendix. 



(67) 



5 Conclusion 

We have considered the one- dimensional problem, assuming that the potential is non- 
negative everywhere. The barrier transmission and reflection amplitudes are described 
in terms of two analytical functions a{K) and /3(k), Eqs. (29-31). Both the functions 
are entire if the potential vanishes outside a finite interval on the x-axis. For potentials 
decreasing exponentially, singularities appear at finite distances to the real axis, corre- 
sponding to the decrease rates at ±00. For «(«;), all the singularities are in the lower 
half-plane, while for the front slope controls the singularities in the upper half- 

plane, and the back slope controls those in the lower half. Poles of the ^-matrix are given 
by zeroes of a(fi;) = 1 — a/2iK. It is proven that for any non-negative potential they are 
all in the lower half-plane. 

The causality in the transmission and reflection processes manifests itself in the an- 
alytical properties of the transition amplitudes. These properties have been employed 
for the space-time description of the tunneling through potential barrier in the Wigner 
phase-space representation p2l . 
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Appendix: Examples 



A number of examples may be found in standard textbooks, e.g. 
i) Square barrier: V{x) = Vq for |x| < Xq, V{x) = for |x| > Xq- 
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sm2a;oP , . 

P = Vo , (68) 

P 

where p — ^/K^ — Vq. One can see that these functions depend, actually, on p^, so there 
is no cut in the K-plane. Zeroes of a are given by (complex) solutions of the equation 

^-^ipxo ^ 1^^^ j . (69) 

This equation has no roots for ImK > 0. If Vq is small enough, there are two roots on the 
imaginary re-axis. Other roots appear in pairs, and their asymptotical position is given 

by 

, TT ^ 1 , 2Rep 
Rep = ±— n. Imp = log —— , (70) 

ZXo Xq Vq 

where n 3> 1 is an integer. 

ii) Exponential barrier: V{x) = Vq exp(— |a;/a;o|). 



[r(i + u)]' ( z 

a — — 



[Jl{z)J_,(z) + Uz)J'_Ml (71) 



sinh 27rre 

where Ji,{z) is the Bessel function, and v — —2iK,, z = 2ixo^/\^. 
iii) The Poschl - Teller barrier: V{x) — Vq/ cosh^(x/xo). 



r2(l - iKXo) ^ _ . COshTTCT ^^^^ 



KXor[-^ + ia — iKXojiy-^ — ia — iKXo) smh ttkxq 

where a — ^^VoXq — |. The function a(«;) has zeroes at kxq — —i{n + |) icr, and (double) 
poles at kxq — —i{n + 1), while 13 has (simple) poles at kxq — ±in, n — 1,2, ■ ■ -. (Note 
that (T is imaginary for 2xo < Vo~^^'^.) 

iv) Narrow barrier: V{x) = Vo6{x). This is the limit one gets as Vq ^ oo, xq — > 0, 
2xoVo = Vq, from the three preceding cases. Now both the entire functions are just 
constant, 

(3 ^vo, (73) 

and a(re) has one zero ai k — —ivo/2. 

v) A double barrier: V{x) = Vi{x — di) +V2{x — ^2)- The case of two non-overlapping 
barriers is described by Eq. (35); a and f3 remain entire functions. As is well known, new 
zeroes may appear in o^k) close to the real axis, corresponding to mctastablc states of 
the particle trapped between the barriers. A special case is that of a symmetrical double 
barrier, where oi = 02 = coshpe"*'' and 61 = —62 = sinhpe*'"'. Now 

a = cosh^ pe~'^'^ + sinh^ pe^*'"^, 
6 = isinh2pcos(K(i -I- 5 + 7), (74) 

where d is the distance between the barrier centers. It is easy to see that the reflection 
may vanish at certain resonance values of the energy, independently of the reflection from 
a single barrier. 
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